The Effect of Shear on Phase-Ordering Dynamics with Order-Parameter-Dependent 

Mobility: The Large-n Limit 
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The effect of shear on the ordering- kinetics of a conserved order-parameter system with 0(n) 
symmetry and order-parameter-dependent mobility T(<t>) ex (1 — <f) 2 /n) a is studied analytically 
within the large-n limit. In the late stage, the structure factor becomes anisotropic and exhibits 
multiscaling behavior with characteristic length scales (£ 2Q+5 /ln£) 1/ ' 2 ( a+2 ' in the flow direction and 
(i/lni) 1 ' 2 '"" 1 " 2 ' in directions perpendicular to the flow. As in the a — case, the structure factor 
in the shear-flow plane has two parallel ridges. 



I. INTRODUCTION 

The dynamics of the ordered phases when a system is 
quenched from the high temperature phase into the low 
temperature region of two- or more-ordered phases has 
been of intense interest W. It is now well established 
that in the late stage, both the equal time pair corre- 
lation function C(r, t) and the structure factor 5(k, t) 
obtained by the Fourier transform of C(r, t) obeys stan- 
dard scaling. By standard scaling it is meant that C(r, t) 
and S(k,t) can be written as f(r/L) and L d g(kL) re- 
spectively, where f(r/L) and g(kL) are the scaling func- 
tions and Lit) is the characteristic length scale in the 
system. For scaling to hold, Lit) must be well separated 
from other length scales that may be present in the sys- 
tem. In most of the systems undergoing phase-ordering 
kinetics, the characteristic length scale Lit) has a power- 
law dependence on the time t elapsed since the quench, 
Lit) ~ t x ' z . The growth exponent z depends on whether 
or not the order parameter is conserved. In the absence 
of shear, but with a constant mobility T, for noncon- 
served order parameter systems z — 2 and for conserved 
order parameter with no hydrodynamic effects z = 3 (for 
scalar fields) and 4 (for vector fields). The character- 
istic length scale L(t) is normally associated with the 
wavevector k m at which the spherically symmetric struc- 
ture factor S(k,t) is maximum (i.e. k m ~ L~~ ). 

When the phase-separating system is subjected to uni- 
form shear (2| , the isotropy in the structure factor is bro- 
ken as the domains grow faster in the flow direction than 
in directions perpendicular to the flow. This anisotropy 
is confirmed by simulation &fl analytical |7J], numerical 
H and experimental results fflll|. We have previously 
shown analytically that within the large-n limit |[7J, the 
structure factor exhibits multiscaling with characteristic 
lengths k~l = k m \ ~ it/hxt) 1 / 4 and k~l ~ it 5 /\nt)^ 
in directions perpendicular to the flow and parallel to the 
flow respectively. For the scalar case (without hydrody- 
namic effects), renormalization group arguments B pre- 
dicts k my ~ t 1 / 3 and k mx ~ i 4 / 3 , for the 'viscous hydro- 
dynamic' regime Lauger et al. |llj found experimentally 

k ml ~ l and Kni ~ * 2 ' Tne ratio kmx/kmy ~ 1/*) IS 

consistent with analytical, numerical, simulational and 
experimental results. 



In the case of unsheared phase-separation, there has 
been much study of the effects of an order-parameter- 
dependent mobility Jl4~|l8t . It was suggested that a mo- 
bility of the form ri<f>) = 1 — 4> 2 is appropriate for deep 
quenches JL2| and to account for the effects of external 
field such as gravity |0|. Simulational calculations E3 ] 
performed on phase-separation with scalar order param- 
eter and mobility r(</>) = 1 — a<j) 2 , show that for a = 1, 
z = A instead of 3 which is a result for constant mobility. 
A crossover from z — 4 to 3 was found for the case where 
a < 1 (with a > 0) [[L4J. For an n- vector order parameter, 
the simulation results Ca] done with r(</>) = 1 — a<p 2 /n 
for n = 2, 3 and 4, shows a crossover from z = 6 to 
z = 4 for a < 1, while for a = 1, z = 6. Emmott et 
al. [ |16|Jl7|| considered a more general expression for the 
mobility given by 



r(0) = r o (i-0 2 /™r, 



(i) 



where a € 5ft + . For scalar fields |L6[ in the Lifshitz- 
Slyozov limit, they found z = 3 + a while for con- 
served vector fields within the large-n limit p7| , mul- 
tiscaling was found with two length scales t 1 / 2 ^' 2 ^ and 
km' 1 ~ (i/ lni ) 1/2(Q+2) - A different form for the mobil- 
ity, T(0) = 1/(1 + exp( a<j) — (3cf) 2 ) ), where both a and (3 
are positive (with (3 > a) was used by Ahluwalia [Q in 
simulating the Cahn-Hilliard model of phase-separation. 
The domain patterns were found to be similar to the ones 
observed in viscoelastic phase separation. It is clear that 
different forms of mobility can be used depending on the 
type of problem concerned. 

In this paper we are mainly concerned with the dy- 
namics of the large-n conserved order parameter with the 
simplest shear flow following a quench from the high tem- 
perature phase to zero temperature, and the case with an 
order-parameter-dependent mobility given by (flh. We 
show that in the late stage, the structure factor becomes 
anisotropic and shows multiscaling behavior ]19| (i.e. 
S(M) ~ iiL x L y L z )^ ( - k - L "' k y L y' k - L ^) with characteris- 
tic length scales L x ~ k mx ~ (i 2Q + 5 /m£) 1/2(Q+2) , L y ~ 

k~l ~ (t/lnt) 1 /2(«+2) ) and Lz „ fc-l „ ( t / lnt )l/2(a+2) 

extracted from the maxima of the structure factor. In 
the k z = plane, there are two parallel ridges, whose 
height and length depends on a. These parallel ridges 



have been observed in experiments M for scalar fields in 
the absence of shear with a = 0. ft is worth mention- 
ing that multiscaling is believed to be an artifact of the 
large-n approximation as in the case of constant mobil- 
ity and zero shear p0[ . For systems with finite n, we 
expect scaling to be recovered asymptotically, although 
multiscaling may be exhibited as a preasymptotic effect 

HI- 

The paper is organised as follows: In the next section, 
model equations which take into the account both shear 
and non constant mobility are introduced. In section III, 
an exact solution for the structure factor is obtained in 
the scaling limit. The discussion of the results for specific 
values of a (i.e a = 1 and 2) are presented in section IV. 
Concluding remarks are given in section V. 



II. MODEL EQUATIONS 

In order to study the phase-separating system, the 
Cahn-Hilliard equation (generalised to n- vector order pa- 
rameter <f>) is given by 



(k x ,k y ,t) — » (k x ,a, t), with t = r and a 



ak x r. 



Or 



r(0)V (- 



V 2 (f> + (j)-{(j) 2 /n)cj) 



(2) 



We are interested in a system with uniform shear flow 
which has a velocity field of the form v = 72/ e x , where 
7 is the constant shear rate and e x is a unit vector in 
the flow direction. For an incompressible system in the 
presence of shear, the term (v • V)0 is added on the left 
hand side of (0) leading to 



dt$ + iyd x (j> 



r(£>v 



-V 2 + 0-(07n)(/> 



In the limit n — ■*■ 00, cj) 2 /n is replaced by its average 
in the usual way, and eq.(||) reduces to a linear self- 
consistent equation whose Fourier transform is given by 

at oky 

where Tq has been absorbed into the time scale, </> is (any) 
one component of <f>, and a(t) = 1 — (<j) 2 )- Eq.(|]) has re- 
cently been solved numerically within a 'self-consistent 
one-loop' approximation for a scalar order parameter 
in two dimensions by Gonnella et al. [p2f. We believe 
that the oscillations (whose amplitudes decreases quite 
rapidly as a increases) found in p2| are slowly-decaying 
preasymptotic transients. 



III. EXACT SOLUTION IN THE SCALING LIMIT 

Eq.(EI) is a first order linear partial differential equa- 
tion which can be easily solved by change of variables: 



With this transformation the left hand side of (|4]) be- 
comes dtfik/dr, and straight forward integration gives 
(after transforming back to original variables) </>k(£) = 
</>k(0) exp f(k, t), where 

/(k, t) = -{n 2 + fi 2 ) 2 b (t) + 4 7 fe x /i (k 2 + f i 2 ) h(t) 
-2 7 2 fc 2 (k 2 + 3^ 2 ) 6a CO + 4 7 3 fc^ ft b 3 (t) 
-7 4 A^ b 4 (t) + (k 2 + ti 2 )Po(t) ~ 2 1 k x ^ Pl (t) 
+J 2 k 2 x P2(t), (5) 

withre 2 = k 2 x +k 2 _, fi = k y +-fk x t, b m (t) = f*dt't' m a(t') a , 

p m {t)=Jldt't' m a{t'Y +l - 

From dimensional analysis, it is easy to see that to 
leading order in t: a(t) ~ t- 1 /(«+ 2 ), L y ~ L z ~ i 1 /2(«+2) 
and L x ~ i(2«+5)/2( Q +2) j with the dominant part of the 
fcz-dependence coming from the shear terms. In fact 
there are logarithmic corrections to these power law rela- 
tions as we will see. It is reasonable to make the ansatz 
a(t) ~ (lni/i) 1 /^ 2 ) (this is true for 7 = [0 ) in the 
large-i limit. Then to leading order in i, we have 



b m (t) = 

Pm(t) = 



2t m bp(t) 

ma + 2m + 2 

t m Po(t) 
ma + 2m + 1 



(6) 



Substituting (1q) in eq. (ra) after making the following 
change of variables, 



jk x 



t 2 b (t) U > Ky 



I pojt) 
bo(t) 



v, k z 



I Po(t) 
bo(t) 



w, 



(7) 



(3) the structure factor S(k,t) = (<^>k(t)^>-k(t)) becomes 



S(M) = Aex-p[2^-F(u,v,w] 

bo 



F(u,v,w) = A (a) -u 3 vAi(a) -2u 2 v 2 A 2 (a) 

5 

o 

-2uv 3 A 3 (a) - v A + —u 2 A 5 (a) 
15 

+— uvAi(a) + v + w — w — 2v w 

o 

2 
- -u 2 w 2 A 2 (a) - 2uvw 2 A 3 (a) , (8) 

o 

where contributions to F which vanish as t — > 00 (at 
fixed u, v, w) have been dropped, A is the size of the ini- 
tial fluctuation, (</>j(r, 0)</>j(r',0)) = A6ij6(r — r'), and 



Ao = 5- 
A x =4- 
A 2 = 'i- 



(>(.) 



40 



40 

a + 4 2a + 6 3a" 
24 24 i 



10 



4a + 10 



a + 4 2a + 6 3a + 8 
12 6 



a + 4 2a + 6 



3 = 2- 


4 


a + 4 ' 


3 


3 


4 2 


2a + 6' 


15 


30 



15 



8a + 24 16a + 40 



(9) 



with A m (a = 0) = 1. In order to find the exact form of bo 
and pq in the large-i limit, we consider the self-consistent 
equation for a(t) given by 



a(t) 



1 



d 3 k 

J2nf 



S(k,t) 



(27r) 3 7 i \b 



Po 



3/2 



x dudv dw exp 2 — F(u, v, w) 
J \ bo 



(10) 



The above integral can easily be evaluated by the method 
of steepest descent using the points of global maxima in 
F(u,V,w) provided Po/bo — * °o as t — ► oo ( this is the 



case we assumed before), therefore eq.(10) becomes 



7tPo 7 V fo o , 



(11) 



where Ci(a) is a constant, F m (o:) is the value at the 
maxima (u m (a),v m (a),w m (a)). The assumption that 
a(t) <C 1 for t 3> 1 has also been used. 



a = 1 


Position (u,v,w) 


Number 


F 


Value 


Type (3D) 


Type (2D) 


(0,0,0) 


1 








Min 


Min 


±(^55/6/3,0,0) 


2 


55/168 


.32738 


IS 


IMax 


±((33 - Vl65)/18V2, -1/V2, 0) 


2 


31/168- ^55/3/28 


.03160 


S2 


S 


±((33 + V165)/18V2, -1/V2, 0) 


2 


31/168+^/55/3/28 


.33744 


SI 


Max 


±(0,0,1/V2) 


2 


1/4 


.25 


SI 


- 


±(V22/3, - v/11/2/4, ±v/5/2/4) 


4 


39/112 


.34821 


Max 


- 


a = 2 


Position (u,v,w) 


Number 


F 


Value 


Type (3D) 


Type (2D) 


(0,0,0) 


1 








Min 


Min 


±(078,0,0) 


2 


14/45 


.31111 


IS 


IMax 


±((7-V7)/4a/2,-1/V2,0) 


2 


31/180 - V7/18 


.02524 


S2 


S 


±((7 + V7)/4V2,-l/V2,0) 


2 


31/180 + V7/18 


.31921 


SI 


Max 


±(0,0,1/V2) 


2 


1/4 


.25 


SI 


- 


±(V35/4,-77/5/2,±V3/5/2) 


4 


59/180 


.32778 


Max 


- 



Table 1. Stationary points of F(u,v,w): Max = maximum, Min = minimum, S = saddle point (2D), Sn = saddle 
point of type n (the matrix of second derivatives has n positive eigenvalues), IS = 'inflection saddle point' (one 
positive, one zero, one negative eigenvalue), IMax = 'inflection maximum' (one zero, one negative eigenvalue). See 
also FIG. 1. for 2D results. 



Therefore to leading logarithmic accuracy, (11) leads 



to 



Pl(t) 



bo In t 3&o hi fro 



2F„ 



2po 



dp db 



dt 



dt 



8F m 

In t 3 In bn 

1 

IF 8F 

^-L m U-i rn. 



(12) 



Using ( p"2| ) and the relations dbo/dt — a", dp /dt = a a+1 
we get 



a(t) 



lnfoo 



1/2 



lni 



-l 1/2 



8F m ln6 32F„ 



(13) 



from which 6p (since dbo/dt = a a ) is found to be 



b (t) 



a + 2 



whence, 



2 lni 

a + 2 



2/(a+2) 



a + 2 



-i N a/ (a+2) 



16F m 32F„ 



pl(t) hit 



b (t) 2F, 
a(t): 



2(a + 2) 



lni 



(a + 2) 2 t 
a + 2 



1/(q+2) 



I V(a+2) 



(14) 



(15) 



16F m 32F, 
The above results for a(t), bo(t) and po(t) justify our orig- 



inal ansatz. From (J7]) we can define the characteristic 
length scales in three directions: L x — 7(£ 2 &o/Po) ~ 
T(a) 7 (t 2Q+5 /lni) 1/2(Q+2) , and L y = L z = (fe /:Po) 1/2 ~ 

T(a)(V lni ) 1/2(Q+2) ; by setting u = fc^x, w = fcy£ v , 
and w — k z L z , where 



T-V) = ^(^ 



2)2 



16F„ 



32F„ 



l/2(a+2) 



(16) 



Exact values for F m (a) and (u m (a),v m (a),w m (a) ) can 
be found by specifying values of a. For example, values 
for a = 1 and 2 are shown in table 1, while for a = 0, 
results presented in [pj are recovered. 

Using (g), (|ll|), and (fl5|), it is easy to show that the 
structure factor becomes 

S(k, t) = Const. [lnU s ] 3/2 Vf (q)/Fm , (17) 

with scaled momentum q and 'scale volume' V s given by 



Q — (fci^n k y L y , k z L z ) , 
U s = L^L 2 ~ 7t (7+2a)/2(a+2) /(lni) 3/2(a+2) 



(18) 



respectively Eq.(|l7|) exhibits multiscaling behavior (i.e 
the power of the 'scale volume' depends continuously on 
the scaling variables). We anticipate, however, that for 
finite n, standard scaling will be recovered and the hit 
terms (which appear in the characteristic length scales) 
will be absent as has been shown explicitly for the case 
with both a and 7 = J20| . 



IV. RESULTS 

Both the positions of the stationary points of F(u, v, w) 
and its value depend on a ( apart from points (0,0,0) and 
±(0, 0, 1/V2) ), see table 1. The global picture of S(k, t) 
is determined by F(u,v,w) (i.e. the stationary points 
and their type together with values of F(u, v, w) at sta- 
tionary points) because In S(k, t) — (F(u, v, w)/F m ) In V s 
[plus k-independent terms] . The features of F(u, v, w) for 
a = 1 and 2 are summarised in table 1 and table 2. 

The parallel ridges found in Figure 1 are similar to 
ones found in experiments |l. The global maxima are 
connected by almost straight ridges to the 'inflexion max- 
ima'. As the value of a increases, the height of each ridge 
decreases towards the limiting value F = 1/4 (and also 
the length of each ridge increases) while the line connect- 
ing the two saddle points and the minimum approaches 
the line u — —v with F = 0. In the (k x ,k y ) plane, the 
ridges become narrower, higher, and closer together as 
a function of k with increasing time t. The angle 6 be- 
tween the ridges and the shear direction (k y direction in 
this case) is a good measure of theory against experi- 
ments because the depth of the temperature quench is 
unimportant as far as the time dependence is concerned. 
We find 



tan(0) = 



C 2 (a) 
It 



(19) 



where C2(a) is a constant which depends on a, e.g. 
C 2 (0) - 2(1-1 A/6), C 2 (2) = (7-V7)/4. Eq. © 
implies that in the (k x , k y ) plane, the ridges move closer 
to the shear direction as time increases, this behavior is 
found both in simulations jq] and experiments fill]. 



F (u, v) 




FIG. 1. The graph of F(u, v, 0) for a — 1. Values for 
F(u,v,0) < —0.1 are not shown. 

In the (u, w) plane, there are four maxima, four sad- 
dle points and the global minimum at the origin. These 
are shown in table 2 (for a — 1,2), and they can eas- 
ily be seen in Figure for a = 1. When a increases, 
the peaks and the higher saddle point (i.e saddle point 
with higher value of F(u, 0, w) reduces towards the lim- 
iting value F = 1/4 (i.e F m (u, 0, w) is deformed towards 
a ring of radius l/\/2)- For each value of a, the structure 
factor pattern in the {k x ,k z ) plane will decrease faster in 
the flow (i.e k x ) direction as a function of t, resulting in an 
elliptical shape with major axis along k z direction. This 
elliptical shape has been observed in experiments [ pT[ for 
scalar fields while the four peaks were not observed. 

The excess viscosity A77, and normal stresses ATVi, 
A7V 2 derived by Onuki || can be evaluated in the asymp- 
totic limit: 

/ffik 
^k x k y S(k,t), 

1 /lnt\ 1/(a+2) 



ry2 W« + 3 

a.Vi = / — (fc 2 -fc 2 )S(M), 



2?r 
In* 



l/(«+2) 



/ii 3 k 
— (k 2 v -k 2 z )S(Kt) 



hit 



l/(a+2) 



(20) 



Numerical calculations [B2| show that both excess vis- 
cosity A-q(t), and the normal stress ANi(t) reach a peak 



before the asymptotic scaling result. It is not possible to 
realise this effect as our calculations are strictly valid in 
the asymptotic regime. 



F (U,W) 




a(t) 



lnA l/(2 +Q ) 



(a + 2) 2 t 

~ a + 2 



1 



16F2P 16F2P 



l/(a+2) 



(22) 



the structure factor can 

2D 



Using both (]2|) and (|22| 

be written down as 5(k, i) = Const. (In j4 s )A s m , 

where A s — L x L y is the 'scale area' and q = 
(k x L x: k y L y ) with L x ~ 7 (i 5+2Q /lni) 1/2(Q+2) , L y ~ 
(t/lnt) 1 / 2 ( Q + 2 ) and F 2D = F(«,u,0). Therefore the 
structure factor pattern is similar to Figure 1. The oscil- 
lations between the peaks G^] which terminate the two 
parallel ridges, we believe are suppressed in the long time 
regime (i.e. they are the preasymptotic decaying tran- 
sients). 



FIG. 2. The graph of F(u,0,w) for a = 1. Values for 
F(u, 0, w) < —0.1 are not shown. 



a = 1 


Position (u,w) 


No. 


F 


Value 


Type 


(0,0) 


1 








Min 


±(^55/6/3,0) 


2 


55/168 


.32738 


S 


±(0,1/V2) 


2 


1/4 


.25 


S 


±(VH0/129,±^/5/43) 


4 


100/301 


0.3322 


Max 


a = 2 


Position (u,w) 


No. 


F 


Value 


Type 


(0,0) 


1 








Min 


±(^7/2/2,0) 


2 


14/45 


.31111 


S 


±(0,1/V2) 


2 


1/4 


.25 


S 


±(^35/3/4, ±1/3) 


4 


17/54 


0.31481 


Max 



Table 2. Stationary points of F(u, 0, w): Max = maxi- 
mum, Min = minimum, S = saddle point. See also FIG. 
2. 

In analogy to 2-D numerical calculations presented in 
p2| , it is easy to repeat the calculation for 2-D in the 
shear-flow (k x , k y ) plane. The self-consistent equation 
for a(t) leads to (with C^,{a) constant) 

1 = -^T exp v ^ m (a) J' (21) 

from which to leading order in t, 



b (t) 



P 2 o(t) 



21nt 
a + 2 



hit 



bo(t) 2F^ 



2/{2+a) 

a + 2 
16F 2D 

1 



1 



\6FZP 



i/(a+2) 



(a + 2) 



V. SUMMARY AND REMARKS 

We have analytically studied the effect of both 
shear and order-parameter-dependent mobility on phase- 
separation within the large-n limit. Shear introduces 
anisotropy in the structure factor pattern because of dif- 
ferent growth rates in the flow direction and directions 
perpendicular to flow. At fixed time t and k, a dis- 
torts the shape of the structure factor S(k, t) (this is 
evident in the F(u,v,0) and F(u,0,w) patterns). Sim- 
ilar to all studies previously done, the order-parameter- 
dependent mobility slows down the rate of coarsening 
(i.e. Li(a — 0) > Li(a ^ 0), where i = x, y, z). 
We believe the multiscaling found here to be the re- 
sult of the large-n approximation, and that for any fi- 
nite n, standard scaling will be obtained with the same 
characteristic length scales but without the lnt terms, 

L x ~ i(2a+5)/2( a +2) ; Ly „ t l/2(a+2) j and ^ „ t l/2(a+2)_ 

The excess viscosity A?y(t), and the normal stresses (i.e 
AJVi(t) and AN 2 (t)) relax to zero as (lni/i Q+3 ) 1/(Q:+2) 
and (lnt/t 1 ^"" 1 " 2 ) respectively in the scaling limit. Again 
we expect logarithmic terms to be absent for finite n. 
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